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Abstract 



We examine the prospects for determining the neutrino mixing matrix 
and for observing CP and T violation in neutrino oscillations in four-neutrino 
models. We focus on a general class of four-neutrino models with two pairs of 
nearly degenerate mass eigenstates separated by approximately 1 eV, which 
can describe the solar, atmospheric and LSND neutrino data. We present a 
general parametrization of these models and discuss in detail the determina- 
tion of the mixing parameters and the mass matrix texture from current and 
future neutrino data in the case where v e and each mix primarily with 
one other neutrino. We find that measurable CP/T violation in long-baseline 
experiments, with amplitude at the level of the LSND signal, is possible given 
current experimental constraints. Also, additional oscillation effects in short- 
and long-baseline experiments may be measurable in many cases. We point 
out that, given separate scales for the mass-squared differences of the solar 
and atmospheric oscillations, observable CP/T violation effects in neutrino 
oscillations signals the existence of a sterile neutrino. We examine several 
textures of the neutrino mass matrix and determine which textures can have 
measurable CP/T violation in neutrino oscillations in long-baseline experi- 
ments. We also briefly discuss some possible origins of the neutrino mass 
terms in straightforward extensions of the Standard Model. 



I. INTRODUCTION 



Our view of the neutrino sector of the Standard Model has recently undergone a rev- 
olutionary change. Observations of solar neutrinos , atmospheric neutrinos , and 
accelerator neutrinos |§ all indicate deviations from their predicted values in the Standard 
Model with massless neutrinos. In each case the observation can be understood in terms of 
neutrino oscillations which in turn requires nondegenerate neutrino masses. The accelerator 
evidence for oscillations is least secure, with preliminary data from the KARMEN Collabo- 



ration JTTJ] excluding some regions of oscillation parameters preferred by the LSND data || . 
Since the solar, atmospheric, and LSND neutrino experiments have different L/E (the ratio 
of oscillation distance to neutrino energy), different orders of magnitude of neutrino mass- 



squared differences 8m are required to properly describe all features of the data [II|. This 
need for three small but distinct mass-squared differences naturally leads to the considera- 
tion of more than three light neutrino flavors. Any additional light neutrino must be sterile, 
i.e., without Standard Model gauge interactions, to be consistent with the well-established 



LEP measurements of Z — > vv HT2|j . From quite general arguments it has been shown |[T3| . |Tl 
that a neutrino spectrum with two pairs of nearly degenerate mass eigenstates, separated 
by a gap of order 1 eV, is required to satisfy all of the constraints from solar, atmospheric, 
accelerator, and reactor data. 

Sterile neutrinos (which we we denote as v s ) have long been considered as an option 
for neutrino oscillations [H)| , |TE|| . More recently a number of models have been proposed 
that utilize one or more sterile neutrinos to describe the existing neutrino data [p] . |r7| -|2~0|| 
or to explain r-process nucleosynthesis pTJ. However, if sterile neutrinos mix with active 
flavor neutrinos they may be stringently constrained by Big Bang nucleosynthesis (BBN). 
In standard BBN phenomenology, the mass-squared difference 5m 2 and the mixing angle 
between a sterile and active neutrino must satisfy the bound 

5m 2 sin 2 29 < 1(T 7 eV 2 , (1) 



to avoid thermal overpopulation of the "extra", sterile neutrino species WM. The restriction 
in Eq. (|l|) would appear to rule out all sterile-active mixing except for small-angle MSW 
or vacuum mixing of solar neutrinos. However, some recent estimates of N u using higher 



inferred abundance of 4 He yield a considerably weaker bound than that given in Eq. ([!]) |23 



Thus BBN may still allow sizeable mixing between sterile and active neutrinos, so models 
with both small and large mixings with sterile neutrinos can be considered. 

In this paper we examine the phenomenological consequences of four-neutrino models in 
which there are two pairs of neutrinos with nearly degenerate mass eigenstates separated 
by about 1 eV, where the mass separations within the pairs are several orders of magnitude 
smaller. We begin with a general parametrization of the four-neutrino mixing matrix, and 
review the current experimental constraints. We then discuss the simple situation where u e 
mixes dominantly with v s or v T in solar neutrino oscillations and mixes dominantly with a 
fourth neutrino [y T or v s ) in atmospheric neutrino oscillations. This situation, which we refer 
to as the dominant mixing scheme, has been shown to fit the existing data reasonably well. 
For dominant mixing we find that the neutrino mixing matrix can be effectively analyzed 
in terms of 2 x 2 blocks, where the diagonal blocks can be approximated by simple two- 
neutrino rotations and the off-diagonal blocks are small but non-vanishing. We then study 
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the relationship between the phenomenology of neutrino oscillations with CP violation and 
the texture of the neutrino mass matrix in models where the two lightest states are much 
lighter than the two heaviest states. 

Since the neutrino mass matrix can in general be complex, and therefore lead to a mix- 
ing matrix with complex elements, CP violation can naturally arise in neutrino oscillations. 
The pertinent question is the size of the violation and how to observe it. We find that if CP 
violation exists, its size may be measurable, and has approximately the same amplitude as 
indicated by the LSND experiment. In many cases there are small amplitude v e — > v T oscil- 
lations that may be measurable in either short- or long-baseline experiments. Furthermore, 
some also have small amplitude —>■ v T oscillations in short-baseline experiments. We dis- 
cuss how oscillation measurements in solar, atmospheric, short- and long-baseline neutrino 
experiments can, in some cases, determine all but one of the four-neutrino mixing matrix 
parameters accessible to oscillation measurements. We also discuss the minimal Higgs bo- 
son spectrum needed to obtain the different types of four-neutrino mass matrices, and their 
consequences for CP violation. 

This paper is organized as follows. In Sec. II we present our parametrization for the 
four-neutrino mixing matrix and expressions for the oscillation probabilities in the case of 
two pairs of nearly degenerate mass eigenstates separated by about 1 eV. We discuss the 
ways in which CP violation, if it exists, may be observed, and the number of observable CP 
violation parameters. In Sec. Ill we summarize the current constraints on the four-neutrino 
mixing matrix and discuss in detail the implications of the dominant mixing scheme. We 
investigate the CP violation effects for several mass matrix textures in Sec. IV. In Sec. V we 
briefly discuss some of the consequences of neutrino mass for non-oscillation experiments, 
such as rare decays and charged lepton electric dipole moments, and we emphasize the 
importance in searching for these rare events, which can reveal new physics effects other 
than neutrino masses. In Sec. VI we summarize our results. Finally, in Appendix A we 
review the number of independent parameters in the mixing matrix of Majorana neutrinos, 
in Appendix B we discuss the modest extensions of the Standard Model Higgs sector that 
allow us to obtain the mass matrix textures, and in Appendix C we determine the neutrino 
mass spectrum and mixing matrix for a particular neutrino mass matrix with CP violation. 



II. OSCILLATION PROBABILITIES 

A. General Formalism 

We work in the basis where the charged lepton mass matrix is diagonal. The most 
general neutrino mass matrix M is Majorana in nature, and may be diagonalized by a 
complex orthogonal transformation into a real diagonal matrix 

M D = U T MU , (2) 

by a unitary matrix U, which is generally obtained from the Hermitian matrix M^M = M*M 
by M^Md = U^M^MU [p4 |. Some general properties of Majorana neutrino mass matrices 



are discussed in more detail in Appendix A. For the four-neutrino case, labeling the flavor 
eigenstates by u x , v e , z/^, v y and the mass eigenstates by vq, V\, v 2 -, ^3, we may write 
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( Vx \ 




fvo\ 




= u 










\Py) 







(3) 



In this paper we will examine the two cases most often considered in the recent literature: 
one of v x and v y is v T and the other is sterile (i/ s ), or both u x and v y are sterile. Explicitly, 
the matrix M may be written in the flavor basis as 



M 



( M xx M xe M Xfl M xy \ 

M xe M ee M efJi M ey 

M Xfl M eM M m 

\M xy M ey M m M yy 



(4) 



The 4x4 unitary matrix U may be parametrized by 6 rotation angles and 6 phases, and 
can be conveniently represented by |25[ 

(5) 



U — R23R13R03R12R02R01 , 



where 



R 



01 



( coi s* 01 QO\ 

-s i coi 
10 
1 / 



(6) 



with 



c jk = cos 6jk 



Sjk = sin 6 jk e 



i5 jk 



(7) 



and the other Rjk are defined similarly for rotations in the j—k plane. The explicit form for 
the 4x4 unitary matrix is 



U 



C01C02C03 


C02C03S01 


~ c 01 0)2^03^13 


— c 02-5oiSo3Sl3 


— CoiCi3So2S^ 2 


— c 13 s 0l s 02Si2 


— c 12Cl3Soi 


+C01C12C13 


— C01C02C13S03S23 


— Co2Cl3So!So3S23 


+CoiSo2Si2 s 13 s 23 


+ s 01 s 02<Si2' S 13 s 23 


— C01C12C23S02 


— Cl2C23 s 01 S 02 


+C12S01S13S23 


— C01C12S13S23 


+C23S01S12 


-C01C23S12 


-C01C02C13C23S03 


-C02C13C23S01S03 


+CQ1C23SQ2SI2S13 


+ c 23 s 01 S 02'Si2'Si3 


+C01C12S02S23 


+ c 12-5oiSo2S23 


+C12C23S01S13 


-C01C12C23S13 


-S01S12S23 


+C01S12S23 



c 03 s 02 



'03 



— s 02 s 03 s 13 c 03 s 13 
+Co2Cl3S 12 



-Ci3S 02 So3S 23 C03C13S23 



+C02C12C23 



— C02C12S23 



(8) 
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We will label the matrix elements of U by U a j , where Greek indices denote flavor eigenstate 
labels (a = x,e,fi,y) and Latin indices denote mass eigenstate labels (j = 0, 1,2,3). With 
the knowledge of the mass eigenvalues and mixing matrix elements, one can invert Eq. 
to obtain the neutrino mass matrix elements 



rrii . 



Pj ut J 



(9) 



3=0 



The vacuum neutrino flavor oscillation probabilities, for an initially produced u a to a 
finally detected up, can be written 



P(v n 



up) = <W ~ E [4 MWlp) sin 2 A kj - 2 hn(W$) sin 2A kj 

j<k 



(10) 



where 



W: 



jk 



a/3 



UajU* k UpjUp k , 



A kj = 5m 2 kj L/(AE) , 5m 2 kj = „u k , 

L is the oscillation distance, and E is the neutrino energy. The quantities W 3 a 
related to the Jarlskog invariants [27] 



2 2 

m h — m A 



(11) 
(12) 

, are 



J aP 



(13) 



which satisfy the identity 



J ap 



jjk* 
J aj3 > 



(14) 



obtained by the interchange U *-> U* in Eq. (|lT|) . Also, 



73 



~ J f3a 



~ J aP ■ 



(15) 



We can also define the real part of Wig as 



Y* = Re(Wi k p) , 



(16) 



which is invariant under interchange of the upper or lower indices: 



Y 



jk 



a 6 



Y, 



kj 



;3a 



Y 



8a 



Y, 



kj 



a0 ■ 



(17) 



Another useful property of the W^p is that the sum over any of the indices reduces them to 
a real positive quantity, e.g., 



ITT I 2 e 
I Uetj | "jk 



a/3 ) 



0. 



(18) 
(19) 



Equations (0), (0), (0D and (0) imply 
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P(u a -> vp) = P{vp -> z/ Q ) , (20) 

which is a statement of CPT invariance. Equation (|T0|) and flTBD imply that nonzero 
can give CP or T violation 

P(z/ Q -> ^) ^ P(z/ Q -> ^) = P(^ -> !/ a ) . (21) 

From Eq. ([H]) we can define the CP- violation quantity 

AP a p = P(v a -> ^) - P(^ -> . (22) 

In four-neutrino oscillations there are only three independent AP a p, and, correspondingly, 
only three of the six phases in U can be measured in neutrino oscillations (for a discus- 
sion, see Appendix A). Thus six angles and three phases can in principle be measured in 
neutrino oscillations, which is the same as in the Dirac neutrino case. Therefore, as far as 
neutrino oscillations are concerned, our results apply equally to Dirac neutrinos. The three 
remaining independent phases in U enter into the mass matrix elements and processes such 
as neutrinoless double beta decay. 



B. Model with two nearly degenerate pairs of neutrinos 

For a four-neutrino model to describe the solar, atmospheric, LSND data and also satisfy 
all other accelerator and reactor limits, it must have two pairs of nearly degenerate mass 
eigenstates ||I3| , |I4|1 ; e.g., Sm 2 un = dm^ <C 5m 2 atm = 5m 2 2 <C ^ m \sND = 8m\ x . We will also 



assume without loss of generality that < m ,mi < m 2 < m 3 . An alternative scenario 
with the roles of 5m% )l and 8m\ 2 reversed gives the same results as far as oscillations are 
concerned, although the implications for the mass matrix, double beta decay and cosmology 
may differ; this alternate possibility will be briefly discussed in Sec. IV. E. Also note that 



if the solar oscillations are driven by the MSW effect [28|, we must require mo > mi; for 
vacuum oscillations, uiq < m\ is also possible. 

Given this hierarchy of the 5m 2 , the oscillation probabilities for a ^ f3 may be written 
approximately as 

P{y a -> vp) ^ ^Isnd sin2 A LS nd + KL sin 2 A atm + B"f m sin 2A atm 

+A^ n sin 2 A sun + B^ n sin 2 A snn , a ^ , (23) 

and for the diagonal channels 

P(u a v a ) ~ 1 - A a L a SND sin 2 A LSND - A™ m sin 2 A atm - A™ n sin 2 A sun , (24) 

where A sca i e = \5m 2 cale L / E , A"{/ ale is the usual CP conserving oscillation amplitude for 
v& — > vp at a given oscillation scale, B^ ale is the CP violation parameter at a given scale, 
and the scale label is sun for the solar neutrino scale, aim for atmospheric and long-baseline 
scales, and LSND for accelerator and short-baseline scales. Note that the CP-violating 
terms have a different dependence on L/E from the CP-conserving terms, which could in 
principle be distinguished by measurements at different L/E ||29|| . In Eqs. (^) and ( p4|) we 
have used the approximation A 31 ~ A 30 ~ A 2 i — A 2 o — A LSND . 
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The oscillation amplitudes are given by 

JaaUp2 + U a3 Up 3 \ 2 = H:\UaQUpQ t u a l u pi 



A&nd = MU a2 UL + U a3 m 3 \ 2 = A\U a0 m + U al U* 01 \ 2 , a?/3, (25) 



Alsnd = KWca? + 


\u a3 


2 )(1- 


-\u a2 \ 2 - 


-\u a3 


2 ), 


= 4(|J7 a0 | 2 + 


\U a x 


2 )(1- 




-\u al 


2 ), 



(26) 

Alt = -4 Re(U a2 U* a3 U*p 2 Up 3 ) , a ^(3, (27) 
A a a ? m = 4\U a2 \ 2 \U a3 \ 2 , (28) 
A*± = -4 Re(U a0 U* al U; U^) , a ? , (29) 
A:: n = 4\U a0 \ 2 \U al \ 2 } (30) 

where the second equality in Eqs. ( p5| ) and (|26|) follows from the unitarity of U. We note that 



the form of the short-baseline oscillation amplitudes in Eqs. ( |25| ) and ( p6|) are different from 
the cases of long-baseline, Eqs. fl2"7]) and (0), and solar, Eqs. ( p9|) and (p0|). The difference 
is due to the fact that the short-baseline oscillations arise from four mass-squared differences 
(5m 20 — 8m\ — ^ m 2 i — ^ m ii)5 while the long-baseline and solar oscillations arise from only 
one mass-squared difference (5m 32 and Sm^, respectively). Probability conservation implies 
Ascaie = ^smie' which can easily be shown using the unitarity of U. The CP violation 

parameters are 

Bit = -2 lm(U a2 U* 3 U*p 2 Up 3 ) , (31) 
Bf un = 2 hn(U a0 U* al U; U pi ) . (32) 

Since B? a = 0, there is no CP violation in diagonal channels. The absence of Bl SND i n 



Eq. ( p3|) shows that no observable CP violation is present for the leading oscillation |30 
and CP violation may only be seen in experiments that probe non- leading scales, Sm 2 tm or 

For short-baseline experiments where only the leading oscillation argument Alsnd has 
had a chance to develop, the off-diagonal vacuum oscillation probabilities are 

P{v a -> vp) ~ A ( f SND sin 2 A L snd , a (3 , (33) 
P(u a -> v a ) ~ 1 - A^jvd sin 2 A L57VD . (34) 

For larger L/.E (such as in atmospheric and long-baseline experiments), where the secondary 
oscillation has had time to develop, the vacuum oscillation probabilities are 

P{y a -> up) ~ ^A^jvz, + A"^ sin 2 A atm + B%f m sin 2A atm , (35) 

P(«/ a - Va ) ~ 1 - - A«f m sin 2 A atm . (36) 

Here we have assumed that the leading oscillation has averaged, i.e., sm 2 A LS nd - ► \- 
Finally, at the solar distance scale, when all oscillation effects have developed, the vacuum 
oscillation probabilities are 

1 



P(u a - up) ~ ^(Af SND + Alt) + At n sin 2 A sun + B^ n sin 2A SW , a ± /3 , (37) 
1 



P{v a - i/ a ) ~ 1 - -(A^ D + A« a J - sin 2 A SWI , (38) 
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where sin 2 A atm has been averaged to | and sin2A atm has been averaged to 0; the CP 
violation at the Sm 2 atm scale is washed out. CP violation is possible only in the off-diagonal 
channels, as noted before, and the solar neutrino v e survival measurement cannot be used 
to observed CP violation. 

More generally, in any model for which the oscillation scales are well-separated and L/E 
is only large enough to probe the largest oscillation scale, CP-violating effects in neutrino 
oscillations will be unobservable (strictly speaking, they are suppressed to order A < 1, 
where A is the oscillation argument for the second-largest oscillation scale) ||30|| . The CP- 
violating effects become observable when L/E is large enough to probe both the largest and 
second-largest oscillation scales. For a four-neutrino model with different oscillation scales 
to describe the solar, atmospheric and LSND data, this means that CP violation can only 
be detected in experiments with L/E at least as large as those found in atmospheric and 
long-baseline experiments. 

As a corollary, in three-neutrino models with two oscillation scales describing only the 
solar and atmospheric data, CP violation has the potential to be observable only in ex- 
periments with L/E comparable to or larger than the solar experiments. However, a mea- 
surement of off-diagonal oscillation probabilities is required to see CP violation, and that is 
not possible in solar neutrino experiments. Hence, if the solar neutrino oscillation scale is 
well-established, the observation of a CP-violation effect in long-baseline experiments could 
imply that there are at least three separate neutrino mass-squared difference scales, and 
thus more than three neutrinos. 



III. DETERMINING THE OSCILLATION PARAMETERS 

In this section we first derive some general constraints imposed by current data on the 
neutrino mixing matrix for four- neutrino models favored by the data, i.e., with two pairs 
of nearly degenerate masses satisfying Sm^ -C Sm 2 2 "C 5m\ v We then determine the form 
of the mixing matrix under the assumption that v e and v x are mostly a mixture of Uq and 
V\ and provide the dominant solar neutrino oscillation, and and v y are mostly a mixture 
of V2 and v 3 and provide the dominant atmospheric neutrino oscillation, which is the form 
of most explicit models in the literature. Then we discuss the measurements needed to 
determine the parameters in the neutrino mixing matrix. The results of this section apply 
equally to the case where rri2 < < mo, mi. The more general case where v e and have 
large mixing with more than one other neutrino is briefly discussed in Sec. IV.E. 



A. Solar v e — > u x and atmospheric z/„ — ► v y 

The flavor eigenstates are related to the mass eigenstates by Eq. Using the formulae 
in Sec. II, the amplitudes for short-baseline oscillation (such as LSND, reactors, and other 
past accelerator oscillation searches) are 

aTsnd = m2u; 2 + u e3 u; 3 \ 2 = A\u e0 u; + u el u;,\ 2 , (39) 
ATsnd = 4(IM 2 + IM 2 )(i - IM 2 - IM 2 ) 

= 4(|f/ M0 | 2 + IM'X 1 - IM 2 - Ki\ 2 ) , (40) 



s 



AT SND =A(\Ue2 


2 + 


|t/e3 


2 )(1- 


l^e 2 | 2 - 


\Ue3\ 


= 4(|[/ e0 


2 + 


\U e i 


2 )(1- 


l^eol 2 - 


Pell 



(41) 



where the second equalities in each case result from the unitarity of U. For atmospheric and 
long-baseline oscillation, the amplitudes are 

A^ m = MU, 2 \ 2 \U, 3 \\ (42) 

AZ m = -4 Re(U e2 U* e3 U; 2 U, 3 ) , (43) 

BZn = "2 lm(U e2 U* e3 U; 2 U, 3 ) , (44) 

AZ m = -4 Re(U e2 U* e3 U; 2 U y3 ) , (45) 

BZn = -2 Im(f/ e2 f/; 3 f/; 2 f/ y 3) , (46) 

A™ m = -4 Re(C/ ft2 C/^C/; 2 C/ y3 ) , (47) 

BZL = -2 \m{U, 2 U; 3 U; 2 U y3 ) , (48) 



and in solar experiments 



AZ n = mo\ 2 \U el \ 2 . (49) 



The atmospheric neutrino experiments favor large mixing of at the atmospheric scale 
Tl,|31fl 



A^>0.8, (50) 
at 90% C.L.; then Eq. ([43J and unitarity imply 

\U l a\ 2 + \U fj3 \ 2 > 0.894, |[/ m0 | 2 + |^i| 2 < 0.106. (51) 
Also, the Bugey reactor constraint [53] gives 

ATsnd < 0-06 , (52) 
over the indicated range for Smj jSND ; then Eq. ( fTi]) implies 

|£4 2 | 2 + |£/ e3 | 2 < 0.016, |t/ e0 | 2 + |t/ e i| 2 > 0.984. (53) 
Finally, the oscillation interpretation of the LSND results || gives 

ATsnd = e 2 , (54) 

where e is experimentally constrained to the range || 

0.05 <e< 0.20, (55) 

where the exact value depends on Sm 2 LSND . Hence, the atmospheric and Bugey results imply 
that \Ue2\i \Uea\, and \U^i\ are all approximately of order e or smaller. We note that 

given these constraints the size of A^ SND must also be small, in agreement with the CDHS 
bound on disappearance |HJ. 
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If we assume that it is only v y that mixes appreciably with in the atmospheric ex- 
periments and only v x that mixes appreciably with v e in the solar experiments, then \U X 2\, 
Uyo\, and |^i| must also be small, i.e., of order e or less. The mixing matrix can 



1^x3 1 j 

therefore be seen to have the form 



U 





u 2 


U 3 


u A 



(56) 



where the Uj (j = 1,2,3,4) are 2x2 matrices and the elements of Ui and Us are at 
most of order e in size. The matrix U4 is approximately the 2x2 maximal mixing matrix 
(i.e., all elements have approximate magnitude l/v2) that describes atmospheric — > v y 
oscillations, but Ui, which is approximately unitary by itself and which primarily describes 
the mixing in the solar neutrino sector, may have large (for vacuum oscillations) or small 
(for MSW oscillations) mixing. 



The form of U in Eq. (|56D, with U2 and U$ ~ e, implies 



(57) 



in the general parametrization of Eq. 
smaller, U takes the form 

/ c i 



After dropping terms second order in e and 



U 



'01 



Q)l 



~ c 0l( s 23 S 03 + C23S02) 
+ s 0l( s 23 s 13 + C23S12) 



" s 0l( S 23 S 03 + C23502) 
" c 0l( s 23 S 13 + c 23 s 12) 



"02 



'12 



C23 



6 03 > 



'13 



"2:-! 



Coi(S23So2 - 
— Soi(s23Si2 



c 23 s 03j 
_ C23S13) 



s 0l(' S 23' S 02 " 
+Coi(s23-Si2 



c 23 s 03j 
" C23S13) 



"S23 c 23 



V 



(58) 



/ 



This matrix provides a general parametrization of the four-neutrino mixing in models where 
u e mixes primarily with v x at the solar mass-squared difference scale, and mixes primarily 
with v y at the atmospheric mass-squared difference scale. Unitarity of U is satisfied to the 
order of e. 

Despite the fact that the expansion of the matrix elements of U in Eq. fl5"8|) is to the first 
order of e, it still allows us, as shown below, to extract all of the interesting oscillation and 
CP violation effects, which are second order in e. Care must be taken when the leading order 
result cancels, and sometimes it is helpful to use the unitarity of U to derive an alternate 
expression that gives the correct leading order answer, e.g., for A^ SND , the first expression 
in Eq. ([25]) gives zero when the form of U in Eq. (|58D is used, but the second expression 
gives a finite (and correct to leading order) result. 

The off-diagonal oscillation amplitudes for the leading oscillation are 



a ex 
LSND 



A efM 

^-LSND 
^LSND 
■rt-LSND 

aw 

LSND 



4|si 2 C 2 3 + S13S23I , 
4|Si2S23 — S13C23I 2 , 
4|s 02 C23 + SosS*23\ 2 , 
0(6 4 ) . 



(59) 
(60) 
(61) 
(62) 
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For atmospheric and long-baseline experiments the oscillation amplitudes are 



Ktm = ~Ktm = -4c 23 Re(st 2 si 3 sy , (63) 

Batm = -Ktm = -2c 23 Im(Si 2 ,Si 3 ,sy , (64) 

A:* m ,B™ m = 0(e 4 ), (65) 

Ktm ^ Ktm = sin 2 29 23 , (66) 

Ktm = -4c 23 Re(so 2 So3sy , (67) 

BaL = 2c 23 Im(so 2 So3sy , (68) 

Batm = -2c2 3 Im[(so 2 so3 + s* 12 s 13 ) s* 23 \ , (69) 



where 823 is defined in Eq. (|7|). All oscillation amplitudes at the LSND scale and all oscillation 
amplitudes (including the CP violation amplitudes) other than A a J m at the atmospheric 
scale are at most of order e 2 . Since A„f m is large, B^ m , which is of order e 2 , may be hard 
to measure since it involves taking the difference of two nearly equal large numbers. At the 
solar scale, we have 

Kin * Kin = sm 2 26 01 . (70) 

There are 6 mixing angles and 3 independent phases that in principle may be measured 
in neutrino oscillations. In all, there are eight independent parameters involved in the 
observables in Eqs. d59|)-([70D, which are the six mixing angles 9 i, 9 02 , 9 03 , # 12 , 6*i 3 , 9 23 , and 
the two phases 

00 = #03 - 5 02 - 5 2 3 , (71) 

01 = £13 - £12 - S 23 . (72) 

These eight parameters could in principle be determined by measurements of the eight ob- 
servables A^ n , Aatm, Aj^ SND , A a t m , B a t m) A^ SND , Aatm, and B^ tm . Therefore if u x = u T , 
then all eight of these parameters could in principle be determined from the solar, atmo- 
spheric, short- and long-baseline experiments. This emphasizes the need for both short- 
and long-baseline measurements of all active oscillation channels, since the oscillation am- 
plitudes involve different combinations of the parameters at short and long baselines. If u x is 
sterile, the three parameters 6*02, #03 an d <fio might be difficult to determine since they would 
involve the disappearance — > u s that is at most of order e 2 in magnitude. If u y = u T , 
the additional observables A e ^ SND , A e a v tmi and B e a \ m can provide a consistency check on the 
parameters # 12 , 6*i 3 , and (f>%. 

We note that from the above results that many of the CP- violating amplitudes can be the 
same order of magnitude as the corresponding CP-conserving amplitudes, and hence poten- 
tially observable in high-statistics long-baseline experiments. The CP violation parameters 
B^tm an d Patm could be determined in vacuum by measuring probability differences AP efl 
and AP UX , where 

AP a/3 = P(y a -> up) - P(u a -> up) , (73) 

in long-baseline experiments, or the probability differences AP efl and AP^ X , where AP a p, 
defined in Eq. (p2|) , measures explicit T- violation. In a vacuum, 
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AP af3 = AP a fs = 2B a a L sin 2A atm . (74) 
Alternatively, one could measure CP asymmetries 

A cp = p K zi A ~ p ^ a ^ A 



or T asymmetries 



a t = P K Zl A ~ P ( U P A (7ft) 
^ P(v a - up) + P{y p - u a ) ' 1 ) 



In vacuum, CPT invariance insures that AP a/3 = AP a p and A^a = A^r '■ However, matter 
effects could induce a nonzero AP a/3 or A^g even in the absence of CP violation P^| . |5T| . 
Since the matter effects in long-baseline experiments for P{y a — > up) and P{yp — > ^ a ) are 
the same, the quantities AP e/1 and ^4^, which can only be nonzero if there is explicit CP 
or T violation, may be preferable |18| . 



There remains a third independent phase that could have consequences for neutrino 
oscillations, but will in practice be difficult to measure. This phase, which could be taken 
as Sqi defined in Eq. @, could be determined from CP violation in u e <-» u^ or u y at the 
5m 2 sun scale, but this effect would require the measurement of an off-diagonal channel at the 
solar scale. Therefore it appears that a complete determination of the four-neutrino mixing 
matrix is not possible with conventional oscillation experiments. Table | lists all parameters 
that appear to be accessible to observation together with the principal observables that 
determine these parameters. 



B. More general mixing scenarios 

In general, both solar u e and atmospheric u^ could oscillate into mixtures of u x and u y . 
In this event #02 and #03 m Eq. (S) are not necessarily small. If one of u x and u y is the 
tau neutrino and the other sterile, there are several possible ways that the existence of such 
mixing could be determined Also, vacuum CP- violation effects involving u e will still be 



no larger than order e 2 (due to the smallness of U e2 and U e3 ), but there are potentially large 
CP-violation effects in long-baseline u^-Uy oscillations (as large as allowed by the unitarity 
of U) p 



IV. CP VIOLATION AND NEUTRINO MASS TEXTURES 



In this section we study the relationship between the neutrino mass texture and the 
possibility for observable CP violation (and, equivalently, T violation) in neutrino oscilla- 
tions in four-neutrino models 
the other is u T (such as in Refs 



We will consider models where one of u x and u y is sterile and 



0, 



and 



and also models where both are 



sterile |37] > which are two possible extensions of the Standard Model neutrinos. Note that in 
all earlier studies the mass matrices were taken to be real and no CP violation was possible. 
In Appendix B we discuss straightforward extensions of the Standard Model for the two 
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cases and show explicitly how their neutrino mass matrices can arise. In Sees. IV.A-IV.D 
we assume that 



m ~ mi <t:m2~m 3 ~^ 5m 2 LSND , (77) 

i.e., the lighter pair of nearly-degenerate mass eigenstates are much lighter than the heavier 
pair, also nearly degenerate, which is the structure of most explicit four-neutrino models in 
the literature. In Sec. IV. E we briefly discuss models with other mass hierarchies. 

From Eqs. (Q) and ( plj| ) the neutrino mass matrix elements involving and v y are, to 
leading order in e, 

M m ~ M* yy ~ m(c 2 2 3 + s 2 23 ) , (78) 
M m ~ im sin 26> 23 sin 5 23 , (79) 
M eil ~ m(s 12 C2 3 + S13S23) , (80) 
M ey ~ m(si 3 c 2 3 - s 12 s* 23 ) , (81) 
M XjJi ~ m(s 2C 23 + s 03 s 23 ) , (82) 
M xy ~ m(s 03 c 23 - s 02 s* 23 ) , (83) 

where we have used the relative sizes of the U a j and mass eigenvalues, and the fact that to 
leading order in e, m 2 ~ m 3 = m. Note that the Sjk, defined in Eq. ([7]), may be complex. 
For the mass matrix elements M xx , M xe , and M ee , all four terms in Eq. @ are small and 
may be of similar size (the first two are suppressed by the small values of mo and mi, the 
last two by mixing angles of size e); their values depend on the exact structure in the solar 
sector, which we do not specify here. More precise solar neutrino measurements would help 
to determine their values. 

The three phases which enter in the expressions for the mass matrix elements given above 
are #23 and 

00 = <5o3 - S 02 + 5 23 , (84) 
4>[ = 813 - 5i2 + 5 23 . (85) 

Only the phases <fto and <pi, which can be measured in oscillation experiments, and <f)' , which 
appears in the expressions for the mass matrix elements, are independent. The two phases 
^23 = (0o ~ 0o)/2 and (/>[ = <p' Q + 4>i — 0o are linearly dependent. 

Equations (|78|)-(|83D may be used to examine the implications of specific textures of the 
neutrino mass matrix. In the following, we discuss several specific textures of the neutrino 
mass matrix which have been considered in the literature. Their CP effects are particularly 
noted. 



A. M, 



eft 



Specific examples of this class of models are given in Refs. [[0J, and ||20|| , in which 
v x — 11 si v y — v ti and the mass matrices are taken to be real. In these models the mass 
matrices were chosen to minimize the number of parameters needed to provide the appro- 
priate phenomenology, and a nonzero M e/Lt is not required. In Ref. [14| the case v x = v T and 



was also considered. 
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Using Eq. (|80|), M eA1 = implies S12C23 — — S13S23, which in turn leads to 

AT SND ~lQ\s 13 \ 2 \s 23 \ 2 sm 2 5 23 , (86) 
AfsND - 4ki 3 | 2 (l - sin 2 26 23 sin 2 5 23 )/c 2 23 , (87) 

Ktm = -AZ m ~ 4|s 13 | 2 |s 2 3rCOs2523, (88) 

BlL = -B e J m ~ -2| Sl3 | 2 |s 23 rsiii2 ( 5 2 3 , (89) 

and 0o = — 25 2 3- Observable oscillations at LSND requires 6*23 7^ 0, 7r and 5 23 7^ 0,n. Also, 
v e —* Vy oscillations in short-baseline experiments, v e — > and v e — > v y oscillations in long- 
baseline experiments, and CP violation in long-baseline experiments are possible, although 
not required, in this scenario. Finally, we have 



|MJ ~ \M yy \ ~ m yjl - sin 2 26 23 sin 2 5 23 , (90) 
\M m \ ~ m sin 26> 23 sin 5 23 . (91) 

As one example, if we take 5 23 — > f and 9 23 — * j, we obtain the model of Ref. |TJL in 
which there is maximal v^-Vy mixing and |M MM | ~ \M yy \ <C |-&^|- Furthermore A e ^ SND = 
Batm — ~Batm — 0, so v e oscillates only to in short-baseline experiments and there is no 
visible CP violation in long-baseline experiments. 



If we allow #23 7^ f , we have a model equivalent to that of Ref. ||19|| ; in this case A^ SND 7^ 



and Batm — ~B e a \ m = 0, so there can be v e — > v y oscillations in short-baseline experiments 
but still no visible CP violation in long-baseline experiments. 

In order to have CP violation in the present case, we must have 523 7^ f • Then there 
must be short-baseline v e — > v y oscillations, although the existence of long-baseline v e <-» z/ M 
oscillations depends on the value of Q 23 . 

Finally, an interesting case to consider is maximal CP violation (maximal in the sense 
that it gives the largest CP- violation parameter for a given |s 13 | and | S23 1 ) , which corresponds 
to #23 = \ ■ If there is also maximal v^-Vy mixing (#23 = j), then the mass matrix in the v^-i/y 
sector is approximately (after appropriate changes of neutrino phase to make the diagonal 
elements real to leading order in e) 

(M m M il y\^m L (li\ 5m 2 atm i7r/4 / 1 1 \ 

[My, M yy )- ^[i l) + Am [l l) ■ [JZ) 
The measurables in short- and long-baseline experiments are then 

ATsnd = A%nd 4kis| 2 , (93) 

A a tm ^atm — j (94) 

t-'atm — t-'atm — l 5 13l > \ JO ) 

i.e., v e oscillates equally into and v y in short-baseline experiments and there are no 
additional contributions to the CP-conserving part of these oscillations in long-baseline 
experiments. The vacuum CP and T asymmetries are especially simple in this case, 

K P = A % = ~ A l - sin 2A atm , (96) 

as the dependence on | S13 1 2 cancels in the ratio. The particular models discussed above are 
summarized in Table |I[ 
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B. M ey = 

An example of this class of models with v x = v s and v y = v T is given in Ref. |36|| , where 
a nonzero M ey was not needed to provide the appropriate phenomenology. From Eq. (j8ip , 
M ey = implies S12S23 — s i3 c 23, which leads to 

ATsnd - 4|s 13 | 2 (l - sin 2 26 23 sin 2 5 23 )/\s 23 \ 2 , (97) 

a lsnd - 16|si 3 | 2 C2 3 sin 2 <5 23 , (98) 

Ktm = -A e Jt m ^ -4ki 3 | 2 c 2 3 cos25 23 , (99) 

B e a L = -B e J m ~ 2| Sl3 | 2 c 2 3 sin25 23 , (100) 

and <fto = —25 23 . The existence of oscillations in LSND implies that 523 7^ f or $23 7^ f • As 
with the M efl = case, it is possible to have v e — > v y oscillations in short-baseline experi- 
ments, v e — > and v e — > v y oscillations in long-baseline experiments, and CP violation in 
long-baseline experiments. The approximate magnitudes of the mass matrix elements M MAt , 
M yy , and M^ y are the same as given in Eqs. (j90|) and (pl|). 



The limit #23 — > and #23 — f reproduces the model in Ref. | 5Efl , which has no 



v P — > v. 



in short-baseline experiments and no visible CP violation in long-baseline experiments. CP 
violation can occur if 523 7^ 0, f ,7T, in which case there are v e — > z/^, oscillations in short- 
baseline experiments and there may be v e — > and v e —>■ v y oscillations in long-baseline 
experiments, depending on the value of 6*23. The M ey = model with maximal CP violation 
and maximal v^-Vy mixing (<5 2 3 = #23 = f ) has the same features as the M eA1 = maximal 
CP violation case in Sec. IV. A, except that the vacuum CP and T asymmetries in Eq. ( |96| ) 
have the opposite sign. The particular M ey = cases discussed here are also summarized in 
Table 0. 



C. M eM + and M ey + 

In this more general case, barring fortuitous cancellations one would expect from 
Eqs. (^[)-(|58]) that there are v e — > and v e —* v y oscillations in short- and long-baseline 
experiments, and CP violation in v e — > and v e — > v y oscillations in long-baseline experi- 
ments. 

For this texture, |M M(U |, \M yy \ <C \M m \ does not necessarily exclude visible CP violation, 
unlike the cases M e „ = or M ey = 0. As an example, the mass matrix 



M 



ei e 2 e l 
e 2 e^ 2 




65 



V e 3 e 



i<P:i 







\ 



e 3 e' 



(101) 



e 6 / 



which is an extension of the model introduced in Ref. [0, leads to CP violation of order e 2 
in v e —>■ and v e —> v y oscillations in long-baseline experiments. This mass matrix differs 
from the one in Ref. 

M MA1 and M yy elements, denoted as 64 and e 6 , respectively, are not necessarily equal, and the 
CP-violating phases are not set to zero. The diagonal elements of the mass matrix can be 



Tl] in that the M efl and elements, denoted as €5, are not zero, the 



15 



taken to be real. Because M ee = M Xfl = M xy = there are only three independent phases. 
The mass eigenvalues and approximate mixing matrix for the mass matrix in Eq. ( |101| ) are 
given in Appendix C. 

The largest off-diagonal short-baseline oscillation amplitudes in this case are 

ATsnd * ±4 , (102) 
AZnd * ±4 • (103) 

Short-baseline v^-Vy oscillations are of order e 4 . The largest long-baseline oscillation proba- 
bilities are 



(104) 

AJtm = ~A a t m ~ e\ — e 3 . (105) 

Short- and long-baseline oscillation amplitudes involving v x are of order e 4 or smaller. The 
observable CP- violating amplitude is 

B e a ? m ~e 3 e 5 sm{(f> 3 + 8 23 ), (106) 

where (see Appendix C) 

, x (e4-e 6 )sin0 1 + 2e 3 e5sin03 

tan<5 23 = — ■ r — — — . 107 

(e 4 + e 6 ) cos <p x + 2e 3 e 5 cos 3 

Finally, 



^ n = sin^ 01 , (108) 



where 



tan #oi 



ei + 2e 3 e 5 ^ ( ei -2e 3 e 5 ) 2 
The phenomenology of these models is summarized in Table O. 



1 | 4eie 3 e 5 (l-cos(0i + 3 -20 2 )) 



D. M ee = M efl = M w = 

In the special class of models with two sterile and two active neutrinos, i.e., both v x 
and v y are sterile, there need not be Majorana mass terms for the two active neutrinos 
in order to obtain the proper phenomenology. As described in more detail in Appendix B, 
M ee = M efJL = M MAt = requires only a minimal extension in the Higgs sector of the Standard 
Model, i.e., only SU(2) singlet Higgs bosons need to be added. Examples of models with 



both v x and v y sterile are given in Ref. |J7]]. We will now show that CP violation effects in 
long-baseline experiments are suppressed in this class of models. 

Models with M eil = have already been discussed in Sec. IV.A; here we add the addi- 
tional constraint = 0. It was previously determined that |M W |, \M yy \ -C \M m \ implied 
5 23 ^ f and 9 23 ^ f , so Eqs. ©-© reduce to 
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A? SND ~ 8\s 13 \ 2 , (110) 
A? SND ^0, (111) 

AZm = -AZm ^ -2| S 1 3 | 2 , (112) 

B% n = -B% n ~0. (113) 

Hence, there is no visible CP violation in long-baseline experiments in this case (strictly 
speaking CP violation is strongly suppressed, to order e 4 ). Therefore in models with M ee = 
M e/1 = M w = 0, the only phenomenological deviations from the Standard Model are CP- 
conserving neutrino oscillations (see Table [H|) an d l ne presence of more than one neutral 
Higgs scalar. 



E. Other mass hierarchies 

There are other mass hierarchies possible which give the same oscillation phenomena as 



those discussed in Sees. IV.A-IV.D. One is m 2 < m 3 <C uiq, m\ ~ \j5m 2 LSND = m, in which 
the solar oscillation occurs between the two upper mass eigenstates and the atmospheric os- 
cillation between the two lower mass eigenstates. Assuming as before that v e mixes primarily 
with v x and with u y (i.e., s 02 , s 03 , s 12 , s 13 ~ e), then 

~ m(c/; [/; + u* el u; x ) • (iu) 

M ee , M ex , M xx , M ey , M Xfl , and M xy are given by similar expressions with appropriate 
changes of subscripts. However, for M MAt , M my and M yy none of the four terms in Eq. @ 
are dominant. Since long-baseline u a — > Up oscillations involve Sm^, and hence the mixing 
matrix elements U a 2, U a3l U^, and (see Eqs. (^2|)-(^)), then a mass texture condition 
such as M efl = 0, when applied to Eq. ( |114| ), does not tell us anything specific about long- 
baseline oscillations. It could, however, affect short-baseline oscillation amplitudes, which 
depend on the U a0 and U a \ (see Eqs. (]3lf)-fl4lD). We do not pursue this possibility further 
here. 

Another possible hierarchy is to have mo — mi < m 2 ~ m 3 where none of the masses are 
much smaller than the others; in this case, all masses would contribute to hot dark matter (an 
alternate possibility, m<i ~ m 3 < m ~ mi with none small, gives similar results). A model 
of this type has been discussed in Ref. pi. In this case, again assuming sq 2 , s 03 , Su, S13 ~ e, 



we find 

M ex ~ m (U: U: + U* el W xl ) , (115) 
with similar expressions for M ee and M xx , and 

M» y ^m 3 (u; 2 u; 2 + u; 3 u; 3 ), (lie) 

with similar expressions for M MAt and M yy . However, for M e/1 , M ey , M Xfl , and M xy , none of 
the four terms in Eq. (|9]) are dominant, the expressions for the mass matrix elements are 
more complicated, and the implications of particular textures for long-baseline oscillations 
are not as easily determined. We also do not pursue this case further here. 
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V. OTHER PROBES OF NEUTRINO MASS 



The presence of mass terms for neutrinos and, in particular, Majorana mass terms, opens 
up a variety of possibilities for phenomena that are not possible in the Standard Model. The 
neutrino mass, besides giving rise to mixing of neutrinos and the associated CP effect dis- 
cussed in this paper, can lead to lepton flavor-changing charged currents analogous to those 
of the quark sector. Majorana mass can also give rise to lepton number violation processes. 
With these possibilities, widely searched-for phenomena such as/i— > e + 7, fx — > e + e + e, 
fi-e conversion, and electric dipole moments for charge leptons, can occur. Unfortunately, all 
these processes [f58[] are proportional to {m u /M w ) A or ((m u / M w )ln{m 2 u / M^)) 41 ■ Given that 
m v is of the order of 5 eV or less p9] these are no larger than 10~ 40 and 10~ 33 . The current 
upper bounds JJDJ are about 30 orders of magnitude larger than those theoretical predictions 



from the neutrino masses. Therefore they are unobservable. Since this conclusion depends 
only on the smallness of the neutrino masses, it is valid in general. 

The Majorana mass term breaks lepton number conservation and can lead to neutrinoless 
double beta decay. The rate is governed by the magnitude of the effective u e mass 

("OHE^iiH^eel, (H7) 
3 

i.e., the magnitude of the M ee element in the Majorana neutrino mass matrix in Eq. (|]). 
The current limit on \M ee \ from neutrinoless double beta decays is about 0.5 eV HTJ]. Since 



M ee = s* lm + (? m m x + sf 2 m 2 + s*^m 3 , (118) 
there will be no visible neutrinoless double beta decay in models with m , mi <C m 2 < m 3 ~ 



5m 2 LSND « 1 eV and \sn\ ~ I s i3| ~ e- In models where m 2 < m 3 <C mo, mi ~ y$ m LSNDi 
or if no neutrino masses are <C 1 eV, neutrinoless double beta decay may provide a strong 
constraint. 



Neutrino masses may also affect cosmology if ^Z v m v > 0.5 eV |[42|| . This level of neu- 
trino mass can easily be accommodated by a four-neutrino model with two pairs of nearly 
degenerate mass eigenstates separated by approximately 1 eV. 

Because of the smallness of the neutrino masses, there are no other observable effects 
besides neutrino oscillations and possibly neutrinoless double beta decay and dark matter. 
However, the rare decays may still be observable if new physics occurs also in other sectors, 
such as anomalous gauge boson couplings or anomalous fermion-gauge boson interactions. 
Therefore, it is important to continue to search for them. 



VI. SUMMARY 

In this paper we have presented a general parametrization of the four-neutrino mixing 
matrix and discussed the oscillation phenomenology for the case of two nearly degenerate 
pairs of mass eigenstates separated from each other by approximately 1 eV, which is the 
mass spectrum indicated by current solar, atmospheric, reactor, and accelerator neutrino 
experiments. We analyzed in detail the case where v e mixes primarily with v x and z/ M with 
u y , where one of u x and u y is v T and the other is sterile, or both are sterile. We found in these 
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cases that the neutrino mixing matrix can be written in 2 x 2 block form with small off- 
diagonal blocks. By construction the mixing matrices have v e — > oscillations in LSND and 
v n — * v y oscillations in atmospheric experiments. We found that the following oscillations are 
also possible: v e — > u y and — > v x in short-baseline experiments, and v e — > i^y,, u e — > i/y f and 
— > v x (including CP- violation effects) in long-baseline experiments. We also found that 
solar, atmospheric, short- and long-baseline oscillation measurements can, in some cases, 
determine all but one of the four-neutrino mixing parameters. Finally, we examined the 
implications of some several specific mass textures, and found the conditions under which 
CP-violation effects are visible. 

As pointed out in the Introduction, additional evidence is needed in order for the four- 
neutrino scenario to be on firm ground. Given that there must be separate mass-squared 
difference scales for the solar and atmospheric oscillations (as currently indicated by the 
data), there are in fact two ways to verify the existence of four light neutrinos: (i) confir- 
mation of the LSND results, which could occur in the future mini-BOONE collaboration 
13|^45[, or (ii) detection of vacuum (i.e., not matter-induced) CP or T violation in long- 



baseline experiments, which should be greatly suppressed in a three-neutrino scenario. 

Once the existence of four neutrinos is established, the next task is to determine the neu- 
trino mixing matrix parameters. We emphasize the significant potential for detecting new 
oscillation channels and CP violation in future high statistics short- and long-baseline oscil- 
lation experiments. Many experiments have been proposed and some will be online in the 
next few years [fi5| , |46|1 . In these experiments neutrino beams are produced at high energy 



accelerators and oscillations can be detected at distant underground detectors. They in- 
clude the KEK-Kamiokande K2K Collaboration IHTfl, the Fermilab-Soudan MINOS [EBfl and 



Emulsion Sandwich |4j| collaborations, and CERN-Gran Sasso ICARUS, Super-ICARUS, 
AQUA-RICH, NICE, NOE and OPERA collaborations Experiments done at muon 

storage rings [0 may be especially important since they will have the ability to measure 
both v e — > and/or z/ T , and — > v e and/or u T , as well as the corresponding oscillation 
channels for antineutrinos. Furthermore, there may also be hitherto undiscovered oscilla- 



tion effects in short-baseline oscillation experiments such as COSMOS |5T[ and TOSCA |p2 
which will search for v T oscillations. To completely determine all accessible parameters 
in the four-neutrino mixing matrix requires searches at both short and long baselines. The 
amplitudes of various oscillation channels, including possible CP violation effects, will help 
further determine the texture of the four-neutrino mass matrix and offer a better under- 
standing of neutrino physics as well as CP violation. 
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APPENDIX A: GENERAL PROPERTIES OF N x N MAJORANA MASS AND 

MIXING MATRICES 



Consider the general case of n neutrinos, in which ur are right-handed (sterile) and ni 
are left-handed (active), % + til = n. We can represent the % right-handed neutrinos by 
their left-handed conjugates and denote the collection of all the n independent left-handed 
neutrinos by a column vector if)^. Then either Dirac or Major ana neutrino mass terms can 
be written as 

^ r M?Ij l + h.c. , (Al) 

with ipn related to ipL by i/jr = ip^C, where C = i^ 2 ^ is the charge conjugation matrix. 
The most general n x n mass matrix M is symmetric and complex, M* = M*, characterized 
by \n{n + 1) magnitudes and the same number of phases. Since a given component field 
in ipL and ipR acquires the same phase factor under a change of phase, n of the phases in 
M may be absorbed into the definitions of the fields, leaving \n(n + 1) magnitudes and 
\n(n — 1) phases; we will often choose the convention that the diagonal elements of M are 
real and the off-diagonal elements complex. 

As shown in Eq. (0), the mass matrix may be diagonalized by a unitary matrix U. 
The effects of U may be divided into two classes of phenomenology which give rise to 
violation of individual lepton number: neutrino oscillations and lepton charged currents. 
Since Z° interacts only with the left-handed neutrinos, neutrino counting in Z — > vv 1 at 
LEP is unchanged if all of the neutrinos are light, i.e., (N U ) LEP = n L , with n L = 3 in 
the Standard Model. This can be demonstrated straightforwardly as follows. The neutral 
current Lagrangian can be written in terms of flavor eigenstates as gipi^KipiZ^ where K 
is an n x n diagonal matrix with the first ur elements being zero and the other elements 
unity. If all neutrinos are very light, which is the case we are considering here, the number 
of neutrinos measured at LEP is simply Tr(KK^) = n^, assuming the couplings of the 
active neutrinos are universal. In terms of the mass eigenstates, the neutrino counting is 
unchanged: Tr(UK(UK)i) = Tr(KK^) = n L . 

In general an n x n unitary matrix such as U can be described by \n(n — 1) rotation 
angles and \n(n + 1) phases. In the charged lepton current, we are free to make phase 
transformations of the charged lepton fields, which removes n of the phases. Then the 
number of surviving measurable phases is \n(n + 1) — n — ^n(n — 1). This argument is 
not affected by the fact that the number of left-handed charged lepton fields, n c , in the 
charged- current is less than n, as each of the first n — n c rows of U can be multiplied by a 
phase without altering the charged currents. Therefore, in general in this Majorana setting 
we can parametrize U by |n(n — 1) angles and |n(n — 1) phases. 

In neutrino oscillations, however, only — l)(n— 2) independent phases can in principle 
be measured, as we will now demonstrate. Note that the W 3 ^ in Eq. (|y]) are invariant when 
U is transformed from either the left or right side by a diagonal matrix which contains only 
phases, i.e., 

U aj -> e^Uaje^ . (A2) 

Then, as far as neutrino oscillations are concerned, we can eliminate 2n — 1 of the phases 
in U, so that there are effectively only 
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-n(n + 1) - (2n - 1) = ~(n - l)(n - 2) , 



(A3) 



independent phases that can be measured by neutrino oscillation experiments. Interestingly 
this is the same number of independent phases that may be determined in the CKM matrix 
for n generations of quarks. However, note that in the most general U there are \n(n — 1) 
phases, so that there are \n(n — 1) — \{n — l)(n — 2) = (n — 1) phases that cannot be 
determined from neutrino oscillations. 

The phase counting of Eq. (|A3|) can also be confirmed by enumerating the number 
of independent CP-violating variables, e.g., the AP a p defined in Eq. (|2"2"|), that can be 
measured. There are \n{n — 1) such differences, but from Eq. ([TJ) AP a/3 = -AP^ and 
from Eq. (|T9| ) J2p AP a p = 0, so it follows that n — 1 of the differences are not independent. 
Therefore there are only |(n — l)(n — 2) independent AP^, and only |(n — l)(n — 2) 
CP-violation parameters can be measured. 



APPENDIX B: HIGGS BOSON ORIGINS OF NEUTRINO MASSES 

The presence of masses for neutrinos is a definite signal of physics beyond the Standard 
Model. Particularly, with the three types of neutrino oscillations which indicate three Sm 2 
scales and require at least four neutrino mass values, a non-trivial extension of the Standard 
Model is necessary. In searching for hints of the extension, it is interesting to consider what 
simplest extensions of the Standard Model are possible and how natural (or unnatural) they 
are in their couplings schemes. In this appendix we discuss briefly the possible origins of 
the two types of mass matrices considered, i.e., models with one or two sterile neutrinos. 
Then both Dirac (active-sterile) and Majorana (active-active or sterile-sterile) neutrino mass 
terms are present. These masses can be obtained by suitable extensions of the Standard 
Model. We will only enlarge the lepton Yukawa sector and the Higgs sector to the extent 
required by the mass matrices of Sec. IV. Our goal is to assure that such mass matrices 
are possible by straightforward modifications of these two sectors, and to determine what 
new particles need to be added to the Standard Model spectrum. We do not attempt to 
construct the best case scenario, which can be done when more information on the neutrinos 
mass are available. For a more extensive discussion of possible origins of neutrino masses 
terms, see Ref. [ f)3|j . 

In the following we denote the right-handed sterile by v S jR an d the corresponding left- 
handed conjugate by P S jL = vJj R C. We also denote the left-handed lepton SU(2) doublet 

by IkL with the corresponding right-handed conjugate fields I^r = —%<J%C\^. R ^ j and k are 
generation labels. 



1. Models with two left-handed and two right-handed neutrinos 

We first consider the class of models in which v x = v s i and v y = u S 2, where v s \ and u S 2 
are the right-handed sterile neutrinos that are associated with v e and u^, respectively, which 
have been considered in Ref. J57]. We discuss two cases: (i) models where only the sterile 



neutrinos have nonzero Majorana mass terms, and (ii) models where both right-handed 
(sterile) and left-handed {y e and v^) fields have Majorana mass terms. In the first case 
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M ee = M efl = M w = 0, which is discussed in Sec. IV. D; in the second case, M e[1 ^ 0, which 
can be realized in the models discussed in Sees. IV. B and IV. C. 

The simplest extension of the Standard Model is case (i) above, which can be obtained 
by adding a singlet real scalar field <fi to the Standard Model Higgs doublet The Majorana 
masses for the two sterile are due to their coupling to <fi and are proportional to the vacuum 
expectation value (vev) of <ft. The $ provides the Dirac masses from Yukawa couplings 
involving both sterile and left-handed neutrinos. We denote the absolute value of the vev's 
of the $ and </> fields as v and V respectively; v is the same as the Standard Model vev. The 
Yukawa couplings can be written as 

2 _ 
£ Y = J2 Gjj'fasjRpsj'L + 9jk{VsjR&lkL + ikR$V sR ) + h.C. , (Bl) 
hi' hk 

where Gjy and g^ are complex couplings, and $ = icr$*. Since u S jR^sj'L = ^sj'L^sjR we 
have Gjji = Gj'j. We also exhibit the symmetry of the Dirac coupling coefficients, g^ = g^j, 

because of the identity v S j R &l kL = IkR^sjR 

For the Higgs potential, we take the simplified case that it has a Z 2 symmetry in <f), i.e., 
invariant under <p — > —0. Then the Higgs potential contains only real coefficients: 

C H = -/i 2 |$| 2 - fi 2 2 (p 2 + A|$| 4 + A 2 4 + A 3 2 |$| 2 . (B2) 

After spontaneous symmetry breaking there are two massive neutral Higgs bosons. Their 
masses are set by v and V. The value of v is the same as in the Standard Model. The 
coefficients of the Yukawa couplings must be very small in order to give neutrino masses 
of the order of eV. In the case V ^> v, one of the Higgs boson is composed mostly of the 
neutral field of $ with a mass proportional to v and the other composed of mostly the </> 
field with a mass proportional to V. There are no other changes to the Standard Model 
phenomenology. 

One can also have a more complicated scenario of case (ii) in which the Majorana masses 
of the two left-handed neutrinos are non- vanishing. These types of models can be constructed 
by the approach discussed below in Appendix B.2. 



2. Models with three left-handed and one right-handed neutrinos 

Here we consider the cases (i) v x = v s and v y = u T , or (ii) v x = v T and u y = u s , where v s 
is a sterile neutrino. In the first case, mass terms are needed in the v ll -v T sector to provide 
the large mixing of atmospheric neutrinos, while in the second case Majorana mass terms 
are needed to provide mixing of solar neutrinos. There are no constraints on which terms 
in the mass matrix may be nonzero, but in each case Majorana masses of the left-handed 
(active) neutrinos must exist. The phenomenology of some of these models is discussed in 
Sees. IV.A-IV.C. 

Majorana mass terms between left-handed neutrinos can arise from the introduction of 



a Higgs triplet which has lepton number —2 [p4| . The Majorana mass of the sterile neutrino 



again comes from a Higgs singlet as discussed above. We denote the triplet by A 
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A= V^ T { *° J ' (B3) 

and the value of the vev of 5° by A. Since A contributes to the masses of the W and Z° 
bosons differently, it has to be small in comparison with the vacuum expectation of the $, 
say A/v < 1CT 2 , so that the bulk of the electroweak gauge boson masses come from the 
doublet. Then, this will not upset the good agreement achieved by the Standard Model 
prediction for the p parameter. 

The Yukawa couplings can be written as 

£y = G s <ph> sR i) sL + QkipsR&lkL + IkL^VsR) + Y h kk ,l kR Al k , L + h.c. , (B4) 

k k,k< 

where G s , g k and h kk i, k and k' = 1,2,3, are complex couplings. The symmetry of the 
Dirac couplings is explicitly exhibited because of the identity v S R&l k L = hh^^sR- h kk i is 
symmetric, h kk i = h k > k , because of the identity l k RAl k iL = l k >RAl k L- For the Higgs potential, 
we can again take the simplified case that is a real scalar field and the Higgs potential has 
the Z 2 symmetry in 0: 

C H = ~ /iil$| 2 - vl<f> 2 ~ ^Tr(AA f ) + r]&A® + r]*&A<<<$> + Ai|$| 4 + A 2 4 

+ A 3 (Tr(AA t )) 2 + A 4 Tr((AA t ) 2 ) + A 5 Tr(A 2 (A t ) 2 ) + A 6 TrA 2 Tr((A t ) 2 ) 

+ «£i0 2 |$| 2 + C20 2 rr(AA f ) + ^ 3 $ t AA t $ + £ 4 &A ] A<5> + £ i5 \$>\ 2 Tr(AA j! ) . (B5) 

The following terms in the Higgs potential, 

r]& A$ + r)*& A f $ , (B6) 

are needed to break the global lepton number invariance in order to avoid the appearance 



of a Goldstone boson called Majoron [p4|| , and r] is the only coupling that potentially can 
be complex. 

To obtain the neutrino mass matrix we can also make phase tranformations on the 
fermion fields v s r and l k i to make G s and the diagonal terms h kk real. If CP is not broken 
spontaneously, which we assume to be the case here, the vacuum expectation values of all 
neutral fields can be made real by phase transformations on the Higgs fields $ and A. After 
spontaneous symmetry breaking this choice of phases for the Yukawa couplings agrees with 
the convention of the neutrino mass matrix discussed in Sec. II. 

With complex couplings, CP violation can generally occur in the Higgs sector. However, 
7] is required to be real by the minimization of the Higgs potential. Hence explicit CP 
violation does not occur in this extended Higgs scenario. A more complicated Higgs potential 
can be chosen to allow complex couplings so that CP violation can be manifest in the Higgs 
sector. We will not elaborate on this possibility here. 

After spontaneous symmetry breaking, the physical Higgs boson spectrum contains a 
doubly charged pair, a singly charged pair, and four neutrals. The Goldstone bosons are 
mostly from the Higgs doublet $. The masses of two of the neutral Higgs bosons are 
proportional to v. The masses of the remaining two neutral Higgs boson are similar to those 
of the case of Appendix B.l. Again the Majorana couplings of the sterile neutrino and the 
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Dirac coupling of the sterile to the left-handed neutrinos are small. However, the Majorana 
couplings among the left-handed leptons do not have to be small if A is chosen to be the 
order of the neutrino masses, i.e., eV ||55|| ; this can be done without leading to any small 
Higgs boson masses. Then the coupling of the singly and doubly charged Higgs boson to 
the charged leptons are not small. The production of these particles in a future high energy 
linear collider or muon collider is possible if they are not too heavy. 

In this type of model, since the constraints M ee = M efl = M MAt = do not apply, the CP- 
violating parameter lm(U e 2U* 3 U* 2 U^) is no longer constrained to be approximately zero. 
The CP violation can be of order e 2 , which is the same order as the oscillation probabilities 
themselves, and hence measurable. Examples of this type of model include the maximal CP 
violation models characterized by Eq. fl92|) and the mass matrix in Eq. ( |101|) . 

We note that to produce the required neutrino masses and the phenomenologically in- 
teresting mass spectrum and expectation values of the Higgs boson fields in both models 
dicussed in this section, new hierarchy problems are introduced |5JJ. In our view, such 
hierarchy problems do not necessarily argue against the models. However, it does argue 
that any model of this sort should be included in a larger, more natural, scheme. Note that 
although the hierarchy in the expectation values of the Higgs boson fields sometimes leads 
to a fine tuning of the parameters, the small vacuum expection value of A may be obtained 
in a natural way p6 . For example, if /X3 3> v 2 , V 2 , rj 2 ^> A 2 , then the minimization of the 
Higgs potential leads to the relation A ~ —r]v 2 /fi 3 <C v. 

There is a growing literature on the generation of neutrino masses. An intriguing class of 
models are those that generate mass dynamically by higher order loop effects |57| . We refer 
the reader to Ref. []5S| for recent and extensive analyses of this possibility for four neutrinos. 
There are also models that use lepton-number violating interactions in i?-parity violating 
supersymmetry for the generation of Majorana mass 15J3 . 



APPENDIX C: AN EXAMPLE WITH CP VIOLATION 

In this appendix we derive the masses and mixing matrix for the model described by the 
mass matrix in Eq. (|101|) . In general a 4 x 4 Majorana mass matrix can have six independent 
phases (see Appendix A), but since three of the mass matrix elements are zero, there are 
only three independent phases in this case. We have chosen to make the diagonal elements 
of M real. 

To achieve the proper neutrino phenomenology, we assume the following hierarchy 

e 2 < ei,e 4 ,e 6 < e 3 ,e 5 < 1. (CI) 
The mass-squared eigenvalues are approximately given by 

2 2 2 2 a 2 2 2 2 /1 1 2 1 2 2\ 2 ( 

m ~ e{m , m 1 ~ 4e 3 e 5 m , m 2 3 ~ (1 + e 3 + e 5 =f e )m , (C2) 

where 

el = 46365 + (e 4 - e 6 ) 2 + 4e 4 e 6 c 2 + 4e 3 e 5 [e 4 cos(0i - <f> 3 ) + e 6 cos(0i + <f> 3 )] , (C3) 

with Cj = cos0j and Sj = sin^-. The eigenvalues are related to the physical mass-squared 
differences by 
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^ m LSJVD 

$ m atm 
,2 
sun 



5m 



2 2 

s r~ 



2 2 I 1 a 2 2\ 2 

m — m 1 ~ (e x — 4e 3 e 5 )m . 



(C4) 
(C5) 
(C6) 



The size of 5m 2 sun in Eq. (|C6l ) implied by the hierarchy of Eq. ( |CT|) means that the solar 
neutrino oscillations are of the MSW type. Therefore in order to have the proper MSW 
enhancement in the sun we must have m\ > mf, which implies |ei| > 1 26365 1. 
The matrix U that diagonalizes M via Eq. (0) is given approximately by 



1001 
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-s 01 e 

essoie^ 01 ^ 3 "* 
V essoie^ 01 ^ 1 ) 



s ie 



-iSoi 



Coi 



-e 5 c ie" 
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-i<t>3 
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V2 



e 5 e 





1 ifea 
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V2 e 



v^ 1 
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where 



tan 5 2 3 
^01 



tan 1 

with c a = cos Q , s a = sin a , a = 03 — 0i 



(e 4 - eg) sin 0i + 2e 3 e 5 g 3 
(e 4 + e 6 ) cos 0i + 2e 3 e 5 c 3 
2e3e5S Q 



ei - 2e 3 e 5 c a 
20 2 , and 



(C7) 

(C8) 
(C9) 



tan 9, 



?2 



01 



1 + 



4eie 3 e 5 (l - c c 



ei + 2e 3 e 5 \ {e l -2e^f 



(CIO) 



We note that this U has the form of of Eq. (|56|). It also can be seen to have the form of 
Eq. (EH!) if we set 



#02 — #03 
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and make the identifications 



sue 



-1013 



V2 
1 

V2 



-i(<^3+<523) 



+ e 5 )e i01 . 



(Cll) 

(C12) 
(C13) 
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TABLES 

TABLE I. Parameters in the four-neutrino mixing matrix and the primary observables used 
to determine them. 
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AflX ,flX rtflX 

LSN Di atm i atm 



TABLE II. Summary of some particular four-neutrino models for rrio,mi <C 777-2 < 777,3 and 
•502) S03) $12, si3 ~ £• All models in the table have been constructed to have short-baseline — > v e 
ocsillations in agreement with the LSND data and large-amplitude — > v y oscillations in atmo- 
spheric and long-baseline experiments; they also all have negligible v e — > u x and — > u y oscillations 
in short-baseline experiments. The size of — > u x oscillations and CP violation in long-baseline 
Vu v y oscillations depend on other model parameters. In all cases, one of v x and v y could be v T 
and the other sterile, or both could be sterile. 
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Sec. IV.D 
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